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ABSTRACT 

The magnetic fields can change their topology through a process known as magnetic reconnection. This 
process in not only important for understanding the origin and evolution of the large-scale magnetic field, but 
is seen as a possibly efficient particle accelerator producing cosmic rays mainly through the first order Fermi 
process. In this work we study the properties of particle acceleration in reconnection zones and show that the 
velocity component parallel to the magnetic field of test particles inserted in nearly non-resistive magnetohy- 
drodynamic (MHD) domains of reconnection without including kinetic effects, such as pressure anisotropy, 
the Hall term, or anomalous effects, increases exponentially. Also, the acceleration of the perpendicular com- 
ponent is always possible in such models. We have found that within contracting magnetic islands or current 
sheets the particles accelerate predominantly through the first order Fermi process, as previously described, 
while outside the current sheets and islands the particles experience mostly drift acceleration due to magnetic 
fields gradients. Considering two dimensional MHD models without a guide field, we find that the parallel 
acceleration stops at some level. This saturation effect is however removed in the presence of an out-of -plane 
guide field or in three dimensional models. Therefore, we stress the importance of the guide field and fully 
three dimensional studies for a complete understanding of the process of particle acceleration in astrophysical 
reconnection environments. 

Subject headings: acceleration of particles — magnetic reconnection — magnetohydrodynamics — methods: 
numerical 



1. INTRODUCTION 

Acceleration of energetic particles is important for a wide 
range of astrophysical environments, from stellar magneto - 
spheres, accretion disk/jet systems, supernova remnants and 
gamma ray bursts to clusters of galaxies. Several mechanisms 
for particle acceleration have been discussed in the literature 
which include varying magnetic fields in compact sources, 
stochastic second order Fermi process in turbulent interstellar 
and intracluster media, and the first order Fermi process be- 
hind shocks. An alternative, less explored mechanism so far, 
involves particle acceleration in magnetic reconnection sites, 
and this will be the focus of the present work which is a first 
of a series of papers on this subject. For a comprehensive re- 
cent revie w on particle acc eleration mechanisms the reader is 
referred to iMelrosel (12009b . 

Magnetic reconnection may occur when two magnetic 
fluxes of opposite polarity encounter each other. In the pres- 
ence of finite magnetic resistivity, the converging magnetic 
lines annihilate at the discontinuity surface and a current sheet 
forms there. It is common knowledge that magnetic fields 
stay frozen in highly conductive fluids. Estimates show that 
Ohmic diffusion for astrophysical scales is absolutely negli- 
gible, contributing to the strongly entrenched view that inter- 
acting magnetic fields of opposite polarity change their topol- 
ogy at low speeds. This is determined by the magnetic re- 
connection speed, and the famous example of two oppositely 
directed magne tic fluxes in Ohmic cont act undergoing slow 
Sweet-Parker (ISweetlll958t iParkerll 19571) process of magnetic 
flux annihilation is frequently invoked. In a more generic con- 
text, the magnetic reconnection rate is the speed at which two 



magnetic flux tubes, which are pushed against each other, can 
pass through altering the initial field topology. This is a situa- 
tion that can be frequently encountered in astrophysical fluids 
with complex motions. 

In the Sweet-Parker model, it has been shown that particles 
can accelerate due to the induc ed electric field in the recon- 
nection zone (Litvinenko 2003). This one -shot acceleration 
process, however, is constrained by the narrow thickness of 
the acceleration zone which has to be larger than the parti- 
cle Larmor radius and by the strength of the magnetic field. 
Therefore, the efficiency of this process is rather limited. Be- 
sides, it also does not predict a power-law spectrum, as gen- 
erally observed for cosmic rays. 

Observations have always been suggestive that magnetic 
reconnection can happen at a high speed in some circum- 
stances, in spite of the theoretical difficulties in explaining 
it. For instance, the phenomenon of solar flares suggests that 
magnetic reconnection should be first slow in order to en- 
sure the accumulation of magnetic flux and then suddenly 
become fast in order to explain the observed fast release of 
energy. A model that can naturally explain this and other ob- 
servation al manifestations of magneti c reconnection was pro- 
posed by Laza rian & Vishniacl (11999b . The model appeals to 
the ubiquitous astrophysical turbulence as a universal trig- 
ger and controller of fast reconnection. The predictions of 
the mod el have been succes sfully tested in numerical simula- 
tions by iKowal et alj (120091) which confirmed that this speed 
is of the order of the Alfven speed in the presence of weakly 
stochastic magnetic field fluctuations. 

An important consequence of fast reconnection of turbu- 
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lent magnetic fields[] is the formation of a thick volume filled 
with reconnected small magnetic flux loops. Now, if such 
turbulent flow is immersed in a current sheet formed by two 
large scale converging magnetic flux tubes (such as, e.g., in 
the Sweet-Parker configuration), then the three-dimensional 
magnetic fluctuations or loops will contract and scatter test 
particles, presenting favorable conditions for energetic parti- 
cle acceleration in a first order Fermi process. In other words, 
the particle may bounce back and forth between these con- 
verging magnetic mirrors formed by oppositely directed mag- 
netic fluxes moving towards each other with the velocity cor- 
respondin g to the reconnection speed. This has bee n first de- 
scribed in Ide Gouveia Dal Pino & Lazarianl (120051) (see also 
lLazariad 12001 " for the situation when there is no back re- 
action of the accel erated partic l es on the reconnected mag- 
netic flu?(Q Later, iDrake etaD (120061) appealed to a similar 
process, but within a collisionless reconnection scenario. In 
their model, the contraction of two-dimensional loops is con- 
trolled by the firehose instability that arises in the particle- 
in-cell (PIC) sim ulations containing both electrons and ions 
(Drak eet al.ll2010l) . In the present work, we show that this 
type of acceleration can be present in a pure magnetohydrody- 
namical scenario as well, where the pressure is fully isotropic 
and the contraction of islands is determined by their interac- 
tions. 

Magnetic reconnection is ubiquitous in astrophysical cir- 
cumstances and therefore it is expected to induce acceleration 
of particles in a wide range of astrophysical environments in- 
cluding galactic and extragalactic ones. For instance, the pro- 
cess has been already disc ussed for the production of ultra- 
high energy cosmic rays (Ide Gouveia Dal Pino & Lazarianl 
120001 120011). a c celera ti on of particle s in ga mma ray bursts 
(Lazarian et al. 2003; Zhan g & Yard 1201 ll) . microquasars 
(Ide Gouveia Dal Pino & Lazarianl 120051) and astrophysical 
jet- acc retion disks in general {de Gouvei aDal Pino et al.l 
l2010alfbl) . In particular, in the case of relativistic jets, a 
diagram of the magnetic energy rate released by violent 
reconnection as a function of the black hole (BH) mass 
spanning 10 9 orders of magnitude was derived and demon- 
strates that the magnetic reconnection power is more than 
sufficient to explain the observed radio outbursts, from mi- 
croquasars to low luminous acti ve galactic nuclei (AGN) 
dde Gouveia Dal Pino et aT1l2010al) . 

More recently, the acceleration in reconnection regions 
has obtained observation al support. It was suggested in 
lLazarian & Opherl (|2009) that anomalous cosmic rays mea- 
sured by Voyagers are, in fact, accelerat ed in the reconnec - 
tion regions of the magnetopause (see also lDrake et al.ll2~010l) . 
Such a model explains why Voyagers did not see any sig- 
natures of acceleration passing the Solar system term i nation 
shock. In a separate development. lLazarian & Desiatil (120101) 
have appealed to the energetic particle acceleration in the 
wake produced as the Solar system moves through interstel- 
lar gas to explain the excess of cosmic rays of the range of 
both sub-TeV and multi-TeV energies in the direction of the 

1 This model of fast reconnection does not impose limits on the amplitude 
of magnetic perturbations. They can be small and the magnetic flux tubes 
require just to have a weak turbulent noise in order to produce fast magnetic 
reconnection, i.e., the reconnection which does not depend on Ohmic resis- 
tivity. 

2 This mechanism is in contrast to the second order Fermi acceleration 
which is frequently discussed in terms of the particle accelerati on by turbu- 
lence generated by reconnection (La Rosa et al. 2006). (see also Kowal et al. 
l20Til in prep.). 



magnetotail. Note, that due to the 1 1 -year solar magnetic cy- 
cle the accumulation of magnetic reversals in the wake is un- 
avoidable. 

The implications of the acceleration process in reconnec- 
tion sites are expected to be even much wider. Numeri- 
cal two-dimension al (2D) simulations recently presented in 
IDrake et ail (120101) confirmed the high efficiency of particle 
acceleration in regions of magnetic reconnection. However, 
we will show here that the process of acceleration happens 
rather differently in two and three dimensional (3D) situa- 
tions. The 3D geometry shows a wider variety of acceleration 
regimes and this calls for much more detailed studies of the 
acceleration. 

In Section [2] we describe the methodology of the particle 
acceleration studies presented here. In Section [3] we show 
the results obtained from these studies, in Sections |4] and [5] 
we discuss the results and draw the main conclusions of this 
work. 

2. METHODOLOGY 

In order to study particle acceleration in magnetic recon- 
nection sites, we performed numerical simulations solving the 
isothermal magnetohydrodynamic (MHD) equations in two 
and three dimensions (2D and 3D, respectively). To compare 
our MHD resul ts with thos e obtained with the Particle in Cell 
(PIC) code in IDrake et alj (120101) we intentionally, in our 2D 
simulations, have reproduced their set up of eight Harris cur- 
rent sheets in a periodic box. The initial density profile is 
chosen in such a way that the total (gas plus magnetic) pres- 
sure is uniform. Initially, we imposed random weak veloc- 
ity fluctuations to this environment in order to enable sponta- 
neous reconnection events and the development of magnetic 
islands. We evolved the system in 2D (both without and with a 
guide field along the third dimension defined by Z direction) 
and in 3D domains. After the development of the magnetic 
loops (which in a 2D geometry form islands) in an underly- 
ing plasma configuration with defined density, velocity and 
magnetic field profiles, we injected test particles at a given 
snapshot and integrated their trajectories solving the equation 
of motion for each charged particle 

d 

— ("/mu) = q(E + uxB), (1) 
dt 

where m, q and u are the particle mass, electric charge and 
velocity, respectively, E and B are the electric and magnetic 

— 1 /2 

fields, respectively, 7 = (l - m 2 /c 2 ) is the Lorentz factor, 
and c is the speed of light. In the studies below we assume 
that the charged particles are protons. 

In the MHD simulations the electric field E is generated 
either by the magnetized plasma or by resistive effects and 
can be obtained directly from the Ohm's law equation, i.e., 

E = -vxB + ry, (2) 

where v is the plasma velocity, j = V x B is the current den- 
sity, and r] is the Ohmic resistivity coefficient. 

In our studies we are not interested in the acceleration by 
the electric field resulting from resistivity effects, therefore we 
have neglected the last term in Equation in the trajectory 
integration. The charged particles then feel an electric field 
v x B besides the magnetic field. Substituting the Ohm's law, 
the equation of motion can be rewritten as 

d 

— (7/wM) = q[(u- v)xB]. (3) 
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Figure 1. Topology of the magnetic field represented as the gray texture with semi-transparent color maps representing locations where the parallel and per- 
pendicular particle velocity components are accelerated for a 2D model with B z = 0.0 at time 6.0 in the code units. The red and green colors correspond to 
regions where either parallel or perpendicular acceleration occurs, respectively, while the yellow color shows locations where both types of acceleration occur. 
The parallel component increases in the contracting islands and in the current sheets as well, while the perpendicular component increases mostly in the regions 
between current sheets. White boxes show regions that are more carefully analyzed in this paper. The simulation was performed with the resolution 8192x4096. 
We injected 10,000 test particles in this snapshot with the initial thermal distribution with a temperature corresponding to the sound speed of the MHD model. 



The particle equation of motion (Eq. [3) was integrated us- 
ing t he 6 th order implicit Runge-Kutta-Gauss (RKG) method 
(see ISanz-Serna & Calvol[T994 e.g.) with a fixed time step 
dt = 10~ 7 . The RKG methods are known to conserve the 
particle energy and momentum in very long integrations 
dSanz-Serna & Cab/olll 994ft in contrast to the standard 4"' or- 
der Runge-Kutta method with the adaptive time step based on 
the 5' h order error estimator (see Pr ess et al.l!992l e.g.), which 
is very commonly used. 

The interpolation of the local values of the plasma velocity 
v and magnetic field B at each step of the integration ha s been 
done using cubic interpolation (Lekien & Marsden 2005) with 
our own discontinuity detector based on a total variation di- 
minishing (TVD) limiter. We performed also the integration 
using linear interpolation which gave us essentially the same 
statistical results, but with much faster integration times. 

In the current study we do not include particle energy 
losses, therefore test particles can gain or lose energy only 
through the interactions with the moving magnetized plasma 
and its fluctuations. The inclusion of radiative and non- 
radiative losses and also the back reaction of the particles on 
the plasma is planned for future studies. 

For simplicity, we assume the speed of light to be 20 times 
the Alfven speed Va, which defines our plasma in a non- 
relativistic regime. The mean density is assumed to be 1 
atomic mass unit per cubic centimeter which is compatible 
with the diffuse interstellar medium (ISM) density. The re- 
sults are presented in units normalized by the assumed light 
speed and time unit (which is 1 hour in our simulations). 

3. RESULTS 



3.1. Parallel and Perpendicular Acceleration in 2D Models 

Figure [ljpresents an evolved 2D configuration of the mag- 
netic field structure with magnetic islands. We clearly see the 
merging of islands in some locations and the resulting defor- 
mations and/or contractions which provide appropriate condi- 
tions for particle acceleration in a similar way to the results 
obtained with the PIC code in dDrake et alJ l2010. and refer- 
ences therein). In general, we distinguish two kinds of accel- 
erating zones where one or both velocity components (parallel 
and perpendicular to the magnetic field) can increase. 

In Figure Q] in order to demonstrate the zones which are 
favorable for a particular type of acceleration, we superim- 
posed on top of the gray texture which represents the mag- 
netic topology, a semitransparent color map showing the lo- 
cations where acceleration increases the parallel (red) or the 
perpendicular (green) velocity component. In yellow areas, 
both types of acceleration can occur depending on the speed 
and direction of the test particle with respect to the orientation 
of the magnetic field lines. 

We see that the increase of the parallel velocity component 
is mostly observed within deformed islands and in current 
sheets (see the red and yellow zones in Figure [TJ, while the 
increase of the perpendicular component is observed mostly 
near and within the islands and between current sheets (see the 
green and yellow zones in Figure [TJ. This complex behavior 
is related to the degree of island deformation and the parti- 
cle direction and speed. If the island contracts, the particle 
can increase its parallel component, however, when the island 
increases its deformation the preferable acceleration usually 
takes place in the perpendicular direction. Interestingly, we 
can find islands in which the acceleration does not happen at 
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Figure 2. Particle energy distributions in the full domain of Figure [T] Blue 
and red histograms show the number of particles accelerating their perpen- 
dicular and parallel velocity components, respectively, at one hour after the 
particles injection in the system. The black line exhibits the initial thermal 
distribution of the injected particles. 

all (see, e.g., the islands above and below the region R3 in 
Figure [TJ. This is due to the fact that these islands are not 
undergoing contraction or deformation processes. 

3.2. Distribution of the Accelerated Particles 

In Figure |2] we present the energy distribution of all parti- 
cles of the domain of Figure Q] at a time corresponding to 1 
hour after the injection, when the particles are still acceler- 
ating. The red and blue lines show the number of particles 
which increase their parallel and perpendicular velocity com- 
ponents, respectively, as a function of the kinetic energy. In 
addition, we show the initial thermal distribution of particles 
(black line). 

In the particle energy distribution corresponding to the per- 
pendicular acceleration component (blue), we distinguish two 
modes, one in the low energy range which is related to the 
initially injected thermal distribution with an extended non- 
thermal tail, and the other corresponding to the high energy 
range. For the particle energy distribution corresponding to 
the parallel acceleration component, we see only one distri- 
bution mode with a maximum at high energies. In the low 
energy range with energies smaller than the particle rest mass 
energy (around 10~ 2 ), there is a clear dominance of the per- 
pendicular acceleration. This creates an anisotropy between 
the two velocity distributions with respect to the magnetic 
field. However, once the particles become relativistic and 
their kinetic energies become comparable to or larger than the 
particle rest mass energy, the anisotropy disappears and the 
acceleration in both directions is equally efficient. 

Figure [3] presents the energy distributions of accelerated 
particles in the selected regions Rl to R6, as shown in Fig- 
ure [T] Regions Rl to R3 correspond to islands where both 
parallel and perpendicular accelerations are observed, while 
regions R4 to R6 correspond to current sheets. 

The energy distributions for regions Rl to R3 are shown 
in the top panels of Figure [3] We notice a dominance of the 
perpendicular acceleration in the low energy range, similarly 
as seen in Figure [2] As we move to higher energies we see 
a trend to form a saturated tail. The efficiency of the accel- 
eration of the parallel component increases with the particle 
energy and at the highest energies the number of particles ac- 
celerating in the parallel and perpendicular directions is com- 



parable. 

In the bottom panels of Figure[3]we present the particle en- 
ergy distributions for regions within current sheets which are 
labeled as R4, R5, and R6 in FigureU In these cases, we also 
observe an anisotropy between the perpendicular and parallel 
components in the low energy range with a dominance of the 
perpendicular acceleration. As in the case of the whole do- 
main, this anisotropy decreases with increasing energy and, 
beyond values around the particle rest mass energy both com- 
ponent distributions become comparable. The high energy 
range tails of these distributions are steeper than those of the 
top diagrams of the figure which correspond to magnetic is- 
lands. This is especially well seen in the energy distribution 
of region R6. This indicates that the number of accumulated 
particles in the high energy range is smaller and therefore the 
acceleration, although still significant, is less efficient in cur- 
rent sheets than within the islands. This is apparently due to 
the fact that (at least in the 2D case) the islands are able to re- 
tain the particles trapped for longer times. What could be the 
dominant mechanism(s) for the acceleration in these zones? 

The acceleration in the contracting/deforming islands is 
suggestive of first order Fermi processes with the parti- 
cles bouncing back and forth between converg i ng m irrors 
as described in Ide Gouveia Dal Pino & Lazarianl (120051) and 
iDrake et alj (120 lOt) . while within and between the current 
sheets the acceleration mechanism is not as clear. In the fol- 
lowing sections, we will examine the acceleration mechanism 
in these different sites in more detail. 

3.3. Acceleration in Contracting Islands 

In Figure |U we show an example of a test proton which is 
trapped in a contracting island and accelerates increasing its 
parallel component. The left panel of the figure shows the 
topology of the magnetic field in a region with a contracting 
island (Rl in Fig. [J) with the trajectory of the proton super- 
imposed. As long as the proton remains trapped in the island 
it orbits around the island center. In the middle panel we plot 
the evolution of the parallel and perpendicular velocity com- 
ponents (red and green lines, respectively), and the kinetic 
energy evolution (blue line) of the test particle. In the right 
panel, we plot the change of the kinetic energy with the par- 
ticle X coordinate. While trapped in the island, the particle 
orbits around the center and increases its energy after each cy- 
cle. Its parallel speed increases while the gyro rotation slows 
down. This results in an exponential growth of the kinetic en- 
ergy of the particle (see middle panel). If we take a closer 
look in the change of the kinetic energy with the position we 
see that the increase of energy happens only when the particle 
moves across the right part of the island. The left panel show- 
ing the magnetic field topology indicates that the field lines in 
this part are contracting due to the interaction with a small is- 
land which is merging with the central one. At the same time, 
the particle gains more energy after each pass in this zone. 

We clearly see how easily the results of Drake e t al.l (1201 Oh 
can be reproduced using the MHD approximation, confirming 
that the process of acceleration in the islands is not restricted 
to the collisionless physics described by PIC codes. MHD 
codes present an easier way to study the physics of particle 
acceleration numerically. 

3.4. Acceleration Near and Within Current Sheets 

We know from shock acceleration theory that particles are 
injected upstream and allowed to convect into the shock, 
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Figure 3. Particle energy distributions in the selected regions of Figure[T] Blue and red histograms correspond to the number of particles which accelerated their 
perpendicular and parallel velocity components, respectively. These histograms are created for all accelerating particles passing through this region from the 
injection moment until 1 hour. 



while diffusing in space so as to effect multiple shock cross- 
ings, and thereby gain energy through the first order Fermi 
process. Now, besides this mechanism, they may also experi- 
ence shock drift acceleration, which is attributed to the grad-B 
drift when the particle encounters an abrupt change in mag- 
netic fields. The origin of this effect is due to the net work 
done on a charge by the Lorentz force (Eq. |2) in a zone of 
non-uniform large scale magnetic field. The principal equa- 
tion governing this is the scalar product of the particle velocity 
(or momentum) and the acceleration by the convective elec- 
tric field, -qv x B. In uniform magnetic fields, the energy 
gain and loss acquired during a gyroperiod exactly cancel, so 
in result no net work is done, AW = 0. 

In contrast, when the gyromotion of a charged particle 
straddles a discontinuity, the sharp field gradient induces an 
asymmetry in the time spent by the charge in either side of the 
discontinuity, so that energy gain and loss do not compensate 
each other. The compressive nature of the field discontinuity 
biases the net work done to positive increments in shock en- 
counters between upstream excursions, and it can be shown 
that AW = qE x dx, where the vxB drifts lie in the x-direction, 
in other words, this energy gain scales linearly with displace- 
ment along the drift coordinate x (i.e., along the discontinu- 
ity). Thus the energy gained by a particle depends on how far 
it drifts along t he front (see, e.g., [Baring & Summerlinl l2009l 
Lugones 2011, and references therein). 

Similarly, in the current sheet zones (regions R4 to R6) 
we may be identifying two distinct mechanisms of accel- 
eration, either first order Fermi acceleration in contract- 
ing/merging islands which are just forming there (or even 
the simple particle scattering between the converging flows 
entering both sides of the current sh eet, as described in 
Ide Gouveia Dal Pino & Lazarianl 12005). or a drift accelera- 
tion, as described above. The islands in current sheets are 
smaller than those of regions Rl, R2 and R3, what results 
in smaller acceleration rates as we see in the lower panels of 
Figure [3] In the zones above and below the current sheets 
it is possible that we see predominantly a drift acceleration 
driven by non-uniformities of the magnetic field. Generally, 



this effect is less efficient than the first order Fermi process 
in merging/contracting islands and results in smaller acceler- 
ation rates. 

In order to better understand these distinct acceleration 
mechanisms, we have also explored the details of the ac- 
celeration of a test particle near and within a single (Sweet- 
Parker shaped) current sheet. Figure |5] shows the trajec- 
tory and energy evolution of this test particle. We see that 
before the particle reaches the current sheet discontinuity 
it is drifted by the plasma inflow and the increasing gradi- 
ent of B as it approaches the current sheet. When it en- 
ters the discontinuity (the white part of the trajectory in the 
left panel), it bounces back and forth several times and gains 
energy (which increases exponentially as shown in the mid- 
dle panel of Figure |5]l due to head-on collisions with the 
converging flow, on both sides of the magnetic disconti- 
nuity (i.e., in a first order Fermi pro cess, as described in 
Ide Gouveia Dal Pino & L azarian 2005). At the same time it 
drifts along the magnetic lines which eventually allow it to 
escape from the acceleration region. Therefore, we see two 
mechanisms: a drift acceleration (dominating outside of the 
current sheet) and first order Fermi acceleration inside the 
current sheet. These processes naturally depend on the initial 
particle gyroradius, since it determines the amount of time the 
particle remains in the acceleration zone before escaping. 

Finally, we may also argue that turbulence above and below 
the current sheets in Figure Q] far from the islands and dif- 
fusion regions, favor second order Fermi acceleration mech- 
anisms with particles being scattered by approaching and 
receding magnetic irregularities. Nevertheless, the first or- 
der Fermi processes occurring within the islands and current 
sheets dominate the overall particle acceleration in the system. 

3.5. The Role of a Guide Field: 2D vs. 3D simulations 

The results presented in the previous sections were obtained 
for 2D models without a guide field. This means that in this 
case the magnetic lines creating the islands are closed and a 
charged particle can be trapped indefinitely in such an island. 
The presence of a guide field normal to the plane of Figure Q] 
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Figure 4. The case of a contracting island where the particles accelerate efficiently (region Rl of Fig. [T}. In the left panel we show the trajectory of a test proton 
trapped in a contracting island. We see two small magnetic islands on both sides of the central elongated island which are merging with it. This process results 
in the contraction of the central island. In the middle panel we show the evolution of the parallel and perpendicular speeds of the test proton (red and green lines, 
respectively) and the evolution of the particle energy (blue line). In the right panel it is shown the change of the particle kinetic energy with the X coordinate. 
The proton orbiting around the center of the magnetic island increases its energy increment after each orbit. 




Figure 5. The case of acceleration near and within a single current sheet with a Sweet-Parker configuration where the particles accelerate efficiently (like region 
R6 of Fig. ff}. The left panel shows the trajectory of a test proton approaching the diffusion region. The color of trajectory corresponds to the particle energy 
(which increases from red to yellow and then finally to white when the particle reaches the current sheet). The middle panel shows the evolution of the parallel 
and perpendicular speeds of the test proton (red and green lines, respectively) and the evolution of the particle energy (blue line). The right panel shows the 
change of the particle kinetic energy with the X coordinate. Outside the current sheet the proton drifts under the effect of the magnetic field gradients and then 
when it arrives in the current sheet it bounces back and forth between the converging magnetic fluxes of opposite polarity while drifting along the magnetic field. 



opens the magnetic loops and allows the charged particles to 
travel freely in the out-of-plane direction. Moreover, the is- 
lands evolve much slower in the presence of a strong guide 
field. 

In Figure© we present the time evolution of the kinetic en- 
ergy of the particles which have their parallel and perpendic- 
ular (red and blue points, respectively) velocity components 
accelerated for three models of reconnection. The upper panel 
shows the energy evolution for a 2D model without the guide 
field (as in the models studied in the previous sections). Ini- 
tially, the particles pre-accelerate by increasing their perpen- 
dicular velocity component only. However, later there is an 
exponential growth of energy mostly due to the acceleration 
of the parallel component which stops after the energy reaches 
values of 10 3 — 10 4 m p (where m p is the proton rest mass en- 
ergy). From that level on, particles accelerate their perpen- 
dicular component only with smaller linear rate in a log-log 
diagram. The middle panel shows the kinetic energy evolu- 
tion in a 2D model with a weak guide field B z =Q.l normal 
to the plane of Figure Q] In this case, there is also an initial 
slow acceleration of the perpendicular component followed 
by the exponential acceleration of the parallel velocity com- 
ponent. However, due to the presence of a weak guide field, 
the parallel component accelerates further to higher energies 
at a similar rate as the perpendicular one. This implies that 
the presence of a guide field removes the restriction seen in 
the 2D model without a guide field and allows the particles 



to increase their parallel velocity components as they travel 
along the guide field, in open loops rather than in confined 
2D islands. This result is reassured by the 3D model in the 
bottom panel of Figure© where no guide field is necessary as 
the MHD domain in fully three-dimensional. In this case, we 
clearly see a continuous increase of both components, which 
suggests that the particle acceleration behavior changes sig- 
nificantly when 3D effects are considered, where open loops 
replace the closed 2D reconnecting islands. 

Considering the parametrization we have chosen for our 
models, the gyroradius of a proton becomes comparable to 
the size of the box when its Lorentz factor reaches a value of 
a few times 10 4 . The largest islands in the system can have 
sizes of a few tenths of the size of the box. These rough esti- 
mates help us to understand the energy evolution in Figure [6] 
and the transition from an exponential to a much slower (lin- 
ear) growth rate in the energy around 10 hours. We note that 
in the diagrams of this figure, the energy is normalized by the 
rest mass value, so that in fact, it is the Lorentz factor that is 
plotted. In the case with absence of a guide field (top panel of 
Fig.©, the exponential parallel acceleration stops right before 
the energy value 10 4 is reached. After this, the rate of accel- 
eration significantly decreases. This occurs because the Lar- 
mor radius of the particles has become larger than the sizes 
of biggest islands. Therefore, from this level on the parti- 
cles cannot be confined anymore within the islands and the 
first order Fermi acceleration ceases. After that, there is a 
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Figure 6. Kinetic energy evolution of a group of 10 protons in 2D models of 
reconnection with a guide field strength B ; =0.0 and 0.1 (top and middle pan- 
els, respectively). In the bottom panel a fully 3D model with initial B-=0.0 is 
presented. The colors show how the parallel (red) and perpendicular (blue) 
components of the particle velocities increase with time. The energy is nor- 
malized by the rest proton mass energy. The background magnetized flow 
with multiple current sheet layers is at time 4.0 in Alfven time units in all 
models. 

much slower drift acceleration (of the perpendicular compo- 
nent only) caused by the gradients of the large scale magnetic 
fields. If a guide field is inserted in such a system, the picture 
is very similar, except for one detail. Now, since the particles 
are able to travel along the guide field, their parallel velocity 
component also continues to increase after the 10 4 threshold 
(see the middle panel of Fig. |6). Of course, in the 3D model, 
the particles follow the same trend (bottom panel of Fig. |6). 

Figure [7] exhibits the dependence of the acceleration rate 
on the out-of-plane guide field strength. It shows the parti- 
cle kinetic energy distribution for three models of reconnec- 
tion with different strengths of the guide field {B z =0.0, 0.5, 
and 1 .0 for the top, middle, and bottom panels, respectively) 
1 hour after the injection. The initial injected thermal distri- 
bution is represented by the blue line. The particles trajectory 
integration was performed at the same MHD snapshot, at time 
4.0 (in units of Alfven time), for all particles in these models. 
We note that at the final state the particle distributions have 
developed high energy tails which depart from the initial ther- 



mal distribution. However, both the number of accelerated 
particles to higher energies and the maximum energy at the fi- 
nal time interval strongly depends on the strength of the guide 
field. 

4. DISCUSSION 

In this paper we have investigated the process of particle 
acceleration in reconnection zones considering test particles 
inserted in an isothermal magnetohydrodynamic domain con- 
taining a series of reconnection layers, without including ki- 
netic effects. There were a few goals that motivated our 
approach. The first goal was to study if the contraction of 
magnetic islands, which develop in 2D reconnection layers 
with finite resistivity, can also occur when the gas pressure is 
isotropic, as is the case in the MH D regime. In earlier stud- 
ies, Drake and collaborators (Dra ke et al.l 120 1 0b invoked ki- 
netic effects such as the firehose and mirror instabilities to 
control the development and contraction of the islands in cur- 
rent sheets. In the present work, after obtaining the formation 
of contracting/deformed islands in current sheets in a nearly 
incompressible MHD domain, the second goal was to inves- 
tigate the acceleration of test particles particularly in the di- 
rection parallel to the large scale magnetic fields during island 
contractions. We were able to identify several regions where 
such acceleration occurs both in a 2D environment with and 
without a guide field and in a 3D environment. Also, we were 
able to identify the nature of the acceleration mechanisms in 
the different acceleration regions, not only within the con- 
tracting islands, but also near and inside the current sheets. 
We found that first order Fermi processes occur in the islands 
and also within the current sheets, while drift acceleration due 
to magnetic field gradients seems to be dominating outside the 
current sheets. 

Our results for particle acceleration i n 2D MHD model s 
were quite similar to those obtained by iDrake et all (120101) . 
i.e., during an island contraction, in our case resulting from 
the merge with other islands, a particle trapped in it can ac- 
celerate and increase its energy exponentially and its parallel 
velocity component grows while the perpendicular one un- 
dergoes a net decrease. During this process the energy gain 
is proportional to the particle energy. These are all necessary 
conditions in a Fermi like process and a close examination of 
the motion of a test particle in a contracting island shows clear 
evidence of the first order Fermi acceleration process. 

Similarly, exponential acceleration of the parallel velocity 
component was also detected in regions within the current 
sheets. A detailed examination of a test particle trajectory 
within and near a current sheet revealed that outside the dis- 
continuity the particle experiences drift acceleration due to 
magnetic field gradients and within the current sheet it goes 
back and forth between the converging flows on both sides 
of the discon tinuity also gaining energy exponentially , just as 
described in Ide Gouveia Dal Pino & Lazarian (2005), while 
d rifting along t he ma gnetic field direction. 

IDrake et alj (|2010) have claimed that particle acceleration 
in an X-point reconnection in a 3D environment, even in the 
presence of a guide field, should behave like in 2D systems. 
However, we have shown that in the presence of an out-of- 
plane guide field , this picture changes completely (see also 
IDrake et al.l2010l) . While in 2D MHD models without a guide 
field the parallel acceleration stops after reaching a certain en- 
ergy threshold, in 2D MHD models with a guide field this 
constraint is removed, and particles can continue increasing 
their parallel speed as they travel along the guide field. Fur- 
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Figure 7. Kinetic energy distribution of the particles for three 2D models 
with different strengths of the out-of-plane guide field B- = 0.0, 0.5, and 1.0 
(top, middle and bottom panels, respectively.) at a time corresponding to 1 
hour after the injection. The initial particle thermal distribution is given by 
the blue line. 



thermore, in fully 3D MHD models with no guide field the 
acceleration of the particles exhibits the same trend as in 2D 
models with a guide field, i.e., there is no constraint on the 
acceleration of the parallel speed. This result is very impor- 
tant as it implies that the overall picture of particle accelera- 
tion in 3D reconnection can be very distinct from that in 2D 
reconnection. Also, it can offer some important clues for par- 
ticle acceleration in more complex domains, such as in the 
presence of 3D turbulence where stochastic reconnection is 



involved (see Kowal et al. 201 1, in prep.). 

Another important result of the present work was to demon- 
strate that the investigated acceleration mechanisms in re- 
connection sites (both drift and Fermi processes) can be 
present in any environment whose evolution can be approx- 
imated by a nearly ideal MHD description as employed here, 
therefore, without the necessity of invoking kinetic instabil- 
ities or anomalous resistivity effects to control the pressure 
anisotropy or the reconnection rate. The nearly non-resistive 
MHD approach with test particle injection offers a possibility 
of exploring more realistic systems. Even though PIC codes 
in general allow the study of plasma processes in greater de- 
tail, they present some disadvantages with regard to an MHD 
description, such as the necessity of a much higher complex- 
ity in the 3D modeling, for i nstance, of turbulenc e. 

A final remark is in order. lOnofri et al.l (120061) investigated 
particle acceleration in reconnection zones and concluded that 
MHD should not be a good approximation to describe the 
whole process of acceleration. However, their 3D numerical 
simulations were performed in a fully resistive MHD regime. 
Therefore, they obtained very efficient particle acceleration 
due to the high electric field induced by resistivity and an 
absorption of most of the available magnetic energy by the 
electrons in a very small fraction of the characteristic time of 
the MHD simulation. This led them to conclude that resistive 
MHD codes are unable to represent the full extent of parti- 
cle acceleration in 3D reconnection. In the present work, we 
have investigated particle acceleration in a nearly ideal MHD 
regime where only numerical resistivity is present. In this 
case, the contribution of a resistivity induced electric field 
is negligible when compared to the advection component, 
namely, the electric field resulting from the plasma motion 
in the magnetized medium, E = -qv x B. Moreover, since 
there is no important dissipation effects, our model is scale 
independent and allows for dimensionless units, so that it is 
straightforward to rescale the configuration to reproduce dif- 
ferent astrophysical environments. The main drawback of the 
numerical resistivity is the creation of a "hole" in the center 
of the magnetic islands, since the numerical modeling cannot 
handle properly the strongly curved magnetic lines over a few 
cells. As a consequence, numerical dissipation removes part 
of the magnetic flux in this region. 

The importance of the present study of acceleration in 
the MHD regime is motivated by the fact that magnetic re- 
connection in 3D becomes fast acco rding to the model i n 
lLazarian & Vishniacl (119991) (see also lLazarian et all 120041) . 
This model has been su c cessf ully tested via numerical sim- 
ulations in iKowal etafl d2009h which confirmed that thick 
reconnection layers form, where magnetic energy is trans- 
ferred into energy of contracting loops 0. In the present pa- 
per, we have confirmed that these loops (or magnetic irreg- 
ularities) can act as the places of efficient particle accelera- 
ti on providing the support for the mecha nism first discussed 
in lde Gouveia Dal Pino & Laz arian (2005). At the same time, 
we see that the process of acceleration does not amount to 
only the process of first order Fermi acceleration outlined in 
the aforementioned paper. Our simulations show a complex 
interplay of different acceleration processes, which motivates 
for further studies of the accele ration in the presen ce of real- 
istic turbulent reconnection (see Kowal et alJ l2009l) . 

The present study has clearly an exploratory character. It 

3 T hese thick layers were also confirmed by Ciaravella & R aymond 
(2008). 
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testifies that the acceleration in reconnection regions may be 
more complex than it may be inferred from earlier simplified 
2D studies. Further study of the acceleration in the MHD 
regime is absolutely essential, as ubiquitous astrophysical tur- 
bulence is expected to induce magnetic reconnection all over 
astrophysical volumes and should be explored in more detail. 
In this situation the acceleration of particles by reconnection 
may play a vital role, the extend o f importance of wh ich can 
be evaluated from further resear ch (Kowal et al.ll20TTh . 

5. SUMMARY 

The results of the paper can be very briefly summarized as 
follows: 

• Advances in the understanding of magnetic recon- 
nection in the MHD reg ime, in particular, tu r bulent 
magnetic reconnection in Lazarian & Vishniac (1999) 
model motivate the studies of whether the reconnection 
in this regime can accelerate energetic particles. 

• Contracting magnetic loops in magnetic reconnection 
in 2D, in the MHD regime, provides the acceleration 
which successfully reproduces the results obtained ear- 
lier with more complicated PIC codes, which proves 
that the acceleration in reconnection regions is a uni- 
versal process which is not determined by the details of 
plasma physics. 

• Acceleration of energetic particles in 2D and 3D shows 
substantial differences, which call for focusing on re- 
alistic 3D geometries of reconnection. Our study also 
shows that apart from the first order Fermi acceleration, 
additional acceleration processes interfere. 
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